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ABSTRACT 

We present the results of a 6.1 square degree survey for clusters of galaxies via their Sunyaev- 
Zel'dovich (SZ) effect at 31 GHz. Prom late 2005 to mid 2007 the Sunyaev-Zel'dovich Array (SZA) 
observed four fields of roughly 1.5 square degrees each. One of the fields shows evidence for significant 
diffuse Galactic emission, and we therefore restrict our analysis to the remaining 4.4 square degrees. 
We estimate the cluster detectability for the survey using mock observations of simulations of clusters 
of galaxies; and determine that, at intermediate redshifts {z ^0.8), the survey is 50% complete to 
a limiting mass (M2oop) of ^ 6.0 x IQ^'^Mq, with the mass limit decreasing at higher redshifts. We 
detect no clusters at a significance greater than 5 times the rms noise level in the maps, and place an 
upper limit on o-g, the amplitude of mass density fluctuations on a scale of 8/i~^ Mpc, of 0.84 + 0.07 
at 95% confidence, where the uncertainty reflects calibration and systematic effects. This result is 
consistent with estimates from other cluster surveys and CMB anisotropy experiments. 
Subject headings: techniques: interferometric, surveys, galaxies: clusters: general; cosmology: observa- 
tions, cosmology: cosmic microwave background, cosmology: large-scale structure 
of the universe 



1. INTRODUCTION 

The number density of massive clusters of galaxies 
depends strongly on cosmology, in particular through 
the normalization of th(^ matter power spectrum, and 
through the dependence of the volume element on the ge- 
ometry of the universe. As the low-energy photons in the 
CMB traverse the hot (~ 10* K) gas of a massive clus- 
ter, about 1% of the photons are inverse-Compton scat- 
tered. The result is a distortion in the CMB spectrum, 
the magnitude of which is proportional to the integrated 
pressure of the intra-clustcr medium (ICM), i.e., the den- 
sity of electrons along the line of the sight, weighted by 
the electron temperature (Sunyaev & Zel'dovich (1972, 
1980); see also Birkinshaw (1999)). The SZ flux of a 
cluster is therefore a measure of its total thermal energy. 

The change in the observed brightness of the CMB due 
to the SZ effect is given by 

= f{x) [ aTUe^^dl = f{x)y (1) 

where Tcmb is the cosmic microwave background tem- 
perature, ctt is the Thomson scattering cross section, ks 
is Boltzmann's constant, c is the speed of light, rrie, ^^e, 
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are the electron mass, number density and temperature, 
respectively, and f{x) contains the frequency dependence 
of the SZ effect (where x = i;^t^^^ ) • Equation 1 defines 
the Compton j/-parameter. The SZ effect appears as a 
temperature decrement at frequencies below « 218 GHz, 
and as an increment at higher frequencies. 

As seen in Equation 1, the ratio of AT/T is indepen- 
dent of the distance to the cluster. This means that the 
SZ effect is redshift-independent in both brightness and 
frequency, offering enormous potential for finding high- 
rcdshift clusters. A cluster catalog resulting from an SZ 
survey of uniform sensitivity has a cluster mass thresh- 
old that is only weakly dependent on redshift for z > 0.7 
(via the angular diameter distance). As a result, SZ clus- 
ter surveys are approximately mass-limited and therefore 
potentially powerful probes of cosmology (e.g., Carlstrom 
et al. 2002). 

Experiments such as the South Pole Telescope (Ruhl 
et al. 2004) and the Atacama Cosmology Telescope 
(Fowler 2004) are surveying hundreds of square degrees of 
sky searching for galaxy clusters through their SZ effect. 
A precursor survey to those being performed by SPT and 
ACT was performed with the Sunyaev-Zel'dovich Array 
(SZA), an 8- telescope interferometer designed specifically 
for detecting the SZ effect towards clusters of galaxies. 
Over the span of two years, the SZA surveyed a small 
region of sky at 31 GHz. This survey has been valu- 
able in characterizing both compact and diffuse cm- wave 
CMB foregrounds. It has resulted in a measurement of 
the power spectrum of the CMB at small scales (Sharp 
et al. 2010), a characterization of extragalactic compact 
source populations (Muchovej et al. 2010), and further 
evidence for large-scale dust-correlated microwave emis- 
sion (Leitch et al. 2010, in prep). In this paper, we 
present results of the survey as they pertain to cosmo- 
logical parameter estimation. 
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The paper is organized as follows: in §2 we describe the 
SZA observations, including a brief description of the in- 
strument, data reduction and calibration, data quality 
tests, and foreground source extraction. We present the 
results of the survey in §3, followed by the calculation 
of the expected number of clusters in §4. In §5 we de- 
termine a constraint on erg and systematic uncertainties 
associated with this analysis are presented in §6. We 
discuss our results in 57. 



2. OBSERVATIONS 



2.1. 



The Sunyaev-Zel'dovich Array 

The Sunyaev-Zel'dovich Array (SZA) is an eight- 
element interferometer located at Caltech's Owens Val- 
ley Radio Observatory near Big Pine, California. The 
observations presented in this analysis were obtained 
with wide-band receivers operating from 27-35 GHz. 
These receivers employ high-electron mobility transis- 
tor (HEMT) amplifiers (Pospieszalski et al. 1995), with 
characteristic receiver temperatures Trx ^ 11 — 20 K. 
Typical system temperatures, including noise due to the 
atmosphere, are on order 40-50 K. Six of the eight anten- 
nas are arranged in a close-packed configuration (spac- 
ings of 4.5-11.5 m), yielding 15 baselines sensitive to ar- 
cminute scales. Two outer antennas yield baselines of 
up to 65 m, for simultaneous detection of contaminating 
compact sources at a resolution of ~ 22.5" (see Muchovej 
et al. (2007) for details of the array configuration). 

This hybrid configuration was chosen to optimize sur- 
veying speed, and can be thought of as a superposition 
of two interferometers: the close-packed array acts as a 
filter for objects of arcminute extent on the sky, such as 
clusters of galaxies, and are used in generating short base- 
line maps; the longer baselines provide high-resolution 
sensitivity to compact sources, and provide the data for 
the long baseline maps. The resulting maps from the 
close-packed array are therefore optimally filtered for 
clusters of galaxies, with radio sources constrained by 
the outlying telescopes. 

2.2. Field Selection 

The SZA survey fields were selected to lie far from 
the plane of the Galaxy and to transit at high elevation 
at the OVRO site, minimizing atmospheric noise while 
optimizing the imaging capabilities of the array. Fields 
were spaced equally in Right Ascension (a) to permit 
continuous observation. These constraints led to the se- 
lection of four regions ranging in declination (6) from 
25° to 35° . No consideration was given to the location of 
known bright radio sources or clusters of galaxies during 
field selection. 

As the SZ effect itself provides no redshift information, 
the fields were also selected to make maximum use of 
publicly available optical survey data, from which photo- 
metric redshifts could be derived. Two of the fields were 
selected to overlap with regions for which optical data is 
available, namely the Deep Lens Survey (DLS Field F2) 
(Wittman et al. 2002) and the NOAO Deep Wide Field 
Survey (NDWFS) (Jannuzi & Dey 1999), also known as 
the Bootes field. Figure 1 depicts the approximate loca- 
tions of the four fields on the IRAS 100 /im dust maps 
(Clegg 1980). 




Fig. 1. — IRAS dust map with overlay of the SZA field locations 



2.3. Data Collection 

Mapping large areas of sky with an interferometer re- 
quires observing multiple discrete pointings within the 
area to be covered, and combining the images from each 
pointing into a mosaic. To perform such a combination 
of images, each of the four fields is split into 16 rows of 16 
pointings. The pointings are equally spaced by 6.6' along 
great circles in the a direction, and each row is equally 
spaced by 2.9' in the S direction. Subsequent rows are 
offset from one another so that that the first pointing in 
each row is shifted by 3.3' in the a direction relative to 
the previous row. This means that for a single field we 
observe an area that spans roughly 2 degrees in the a 
direction and 1 degree in the 6 direction (see Figure 2). 




Fig. 2. — Mosaic pointing locations for a given SZA survey field. 
The fields are divided into 16 rows of 16 columns, with the pointings 
in each row separated by 6.6' and each row offset from each other 
by 2.9'. This leads to each field being roughly 2 degrees by 1 degree 
in area. In a single track the SZA observed four pointings within a 
given row. For example, pointings in the first and ninth, followed 
by pointings in the second and tenth columns. 



For each of the survey fields, data were taken daily 
in 6 hour tracks. In a single track, we observed two 
staggered pairs of pointings, within a single row. These 
observations were performed in a manner that permits 
ground subtraction from consecutive pointings in a pair 
(although the ground contamination was found to be neg- 
ligible and was not subtracted as part of the analysis of 
this paper) . Each track results in roughly 1 hour of ob- 
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servation on each of the four pointings, with very nearly 
the same Fourier sampling for pairs of pointings. A sec- 
ond track is run at a later date, with the order of the 
pairs reversed, to ensure that the Fourier sampling for 
all four pointings is comparable. In Figure 2 we show 
the position of the pointings in each field, and indicate 
how the pointings were observed in a given track. 

For each set of four pointings, this sequence is repeated 
three times over the span of roughly one year, so that 
each pointing is observed in six tracks, translating to 
roughly 6 hours of observation per pointing over the du- 
ration of the survey. 

2.4. Observations and Data Reduction 

Images of the survey fields were produced by linear 
mosaicking of maps from the individual pointings, as de- 
scribed in a companion paper, Muchovej et al. (2010), 
hereafter MIO. In particular, we stitch together maps 
from the individual pointings, properly weighted by the 
primary beam, to generate signal and noise maps of the 
fields. We further construct a significance {snr) map by 
taking the ratio of the signal and noise maps. In Ta- 
ble 1 we present details of the mosaicked SZA survey. 
The second and third columns show the approximate 
center of each 16-row field. Wc also present the band- 
pass and gain calibrators in the next two columns, with 
their fluxes as measured by the SZA (calibrated to ob- 
servations of Mars, see below). In the fifth column we 
give the time range over which observations were taken, 
with the caveat that observations were not performed 
every day during that time span. The penultimate col- 
umn lists the total integration time for data used in the 
analysis, and the final column gives the number of rows 
observed in each field. To ensure uniform coverage of all 
fields, tracks were repeated when deemed necessary due 
to poor weather or instrumental glitches. Note that the 
full 16 rows were not observed for all fields, due to main- 
tenance operations, instrumental characterization, and 
RFI monitoring. For the first 8 months of observations, 
the SZA4 field was dedicated to CMB anisotropy mea- 
surements (Sharp et al. 2010), and data were collected 
in a manner incompatible with survey observations. As 
a result, only 7 rows in the SZA4 field were completed to 
the full survey depth. We observed 6 more rows in that 
field, but not to the same depth as the rest of the survey. 
This results in a smaller region of uniform sensitivity in 
SZA4, but the field is still usable in our analysis. In 
total, the data in the SZA cluster survey correspond to 
1493 tracks taken between November 13, 2005 and July 
25, 2007. 

Data for each track were calibrated using a suite 
of MATLAB^ routines, which constitute a complete 
pipeline for flagging, calibrating, and reducing visibil- 
ity data (see Muchovej et al. 2007). Although the data 
were reduced exactly as described in that work, survey 
data collection differed in a few key ways: whereas in 
targcited observations we observed a source for 15 min- 
utes before observing a calibrator, in this work four dis- 
tinct pointings were observed for roughly 4 minutes each 
before observing a calibrator. Also, system temperature 
measurements were performed every eight minutes in sur- 

9 The Mathworks, Version 7.0.4 (R14), 

http : //www . mathworks . com/products/matlab 



vey mode, as opposed to every 15 minutes in targeted 
observations. The absolute flux calibration is referenced 
to Mars, assuming the Rudy (1987) temperature model. 
Accounting for the uncertainty in the Mars model and 
in the transfer to our data, we assign a conservative un- 
certainty of 10% to our flux calibration. Flagging of the 
data as described in Muchovej et al. (2007) resulted in a 
loss of roughly 23%. At the end of a single 6-hour track, 
our on-source time per pointing was roughly 55 minutes, 
leading to a noise level of approximately 1.5 mJy/beam 
in each pointing of the short and long baseline maps. 
Lastly we remove approximately 0.4% of the data with 
poor noise properties resulting from minor glitches in the 
digital correlator. 

2.5. Mosaics 

Once data on all pointings in a given field are reduced, 

we construct a linear mosaic of the field on a regular 
grid of 3.3" resolution. This scale is much less than the 
requirement for Nyquist sampling of the data, lyfY" — ) 
where D^ax is the longest baseline, and leads to a con- 
venient number of pixels for the use of FFTs in the fol- 
lowing analysis. The maps arc composed of the data 
across our 8 GHz of bandwidth, centered at a frequency 
of 30.938 GHz. The attenuation due to the primary 
beam response for each pointing is corrected before mo- 
saics are constructed. The primary beam is calculated 
from the Fourier transform of the aperture illumination 
of each telescope at the central observing frequency, mod- 
eled as a truncated Gaussian with a central obscuration 
corresponding to the secondary mirror. Typical synthe- 
sized beams for the short and long baseline maps for 
each of our pointings have Gaussian FWHM of 2' and 
45" FWHM, respectively. 

Due to the overlap of neighboring pointings, the ef- 
fective noise is approximately uniform in the interior of 
the mosaics, but increases significantly towards the edge 
of the mosaicked images. Wc limit the survey area by 
applying an edge cutofl' in our mosaicked maps where 
the; effcictive noise is > 0.75 rnjy/beam (corresponding 
roughly to the one-third power point of the beam, given 
the noise in a single pointing). 

In Table 2 we show the noise properties of the observed 
fields. Wc present the minimum and median noise (in 
mJy/beam) for mosaic maps made with long baselines 
only, short baselines only, and with the combination of 
the two. The median noise is calculated only in the region 
within which the noise is less than the 0.75 mJy/beam 
cutoff. The last column indicates the total area covered 
in each field below the noise threshold. That the min- 
imum and median pixel noise values are similar is an 
indication of the uniformity of the coverage in the sur- 
vey fields. The SZA4 field is not as uniform as the other 
fields, as only seven rows were completed to full survey 
depth, while the remaining six rows were observed for 
roughly half the time. 

2.6. Source Extraction 

Extra-galactic radio sources are a significant contami- 
nant in observations of the SZ effect, particularly at cen- 
timeter wavelengths. The source extraction algorithm for 
the SZA survey consists of two stages, both of which use 
5 GHz VLA follow-up observations of our fields to assist 
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TABLE 1 
Survey Observations 



Field Name 


Field Center 


(J2000) 


Calibrators 


Dates 


Integration 


Rows 




a 


5 


Bandpass (Jy)^ 


Gain (Jy)» 


of Observations 


Time (hrs) 


Covered 


SZA4 


02''15'"38''.3 


32°08'21" 


J2253+161 (11.6) 


J0237+288 (2.9) 


07/11/2006 to 07/25/2007 


687 


7 


DLS 


09'' 19"* 40°. 


30°01'26" 


J0319+415 (11.0) 


J0854+201 (5.4) 


11/18/2005 to 07/06/2007 


1054 


14 


NDWFS 


14'' 30"* 08°. 


35°08'34" 


J1229+020 (25.3) 


J1331+305 (2.1) 


11/19/2005 to 07/23/2007 


1000 


14 


SZA3 


2l''30"'07°.0 


25°01'26" 


J1642+398 ( 5.5) 


J2139+143 (1.4) 


11/13/2005 to 07/25/2007 


1245 


16 



^Fluxes obtained from 31 GHz SZA observations of sources on April 16, 2006. 



TABLE 2 
Survey Sensitivity 



Field 


Short Baselines 


Long B 


laselines 


AH Baselines 


Area 


Name 


Minimum rms 


Median rms 


Minimum rms 


Median rms 


Minimum rms 


Median rms 


Covered 




(mjy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(degree^) 


SZA4 


0.218 


0.400 


0.231 


0.422 


0.159 


0.305 


1.5 


DLS 


0.201 


0.237 


0.200 


0.250 


0.142 


0.173 


1.5 


NDWFS 


0.219 


0.239 


0.219 


0.247 


0.156 


0.172 


1.5 


SZA3 


0.213 


0.232 


0.218 


0.241 


0.153 


0.167 


1.7 



in source identification. The first stage is an iterative fit- 
ting of bright sources with fluxes at least 5 times the local 
map rms. This stage of fitting is described extensively in 
MIO and summarized in this section. The second stage 
of source removal relies heavily on VLA follow-up obser- 
vations to remove dim sources. A description of the VLA 
follow-up observations, and their use in determining pa- 
rameters of the sources to be fit in the SZA data, can 
also be found in MIO. 

Using the following algorithm we fit for 326 total 
sources: 239 in the iterative stage of fitting (see §2.6.1), 
and 87 in the second stage of fitting (see §2.6.2). We 
note that of the 326 sources, 39 (12%) were deemed to 
be extended in the follow-up VLA data and fit as ex- 
tended sources. However, the fitted 31 GHz major axes 
were all determined to be smaller than 22.5" (the FWHM 
of the long baseline maps), implying they are mostly un- 
resolved by the SZA. This allows us to approximate all 
sources detected in our survey as unresolved by our long 
baselines in subsequent analysis (see §4.2). 

2.6.1. First Stage: Bright source removal 

The first stage of the source extraction algorithm con- 
sists of iterative removal of bright sources. Source iden- 
tification begins in the image plane, with inspection of 
the combined (short and long baseline) significance (snr) 
maps for the brightest pixel with significance greater 
than 5. Once we identify the location of a source, we next 
determine whether the source is extended or unresolved 
as seen by the VLA, and whether this candidate is a sin- 
gle source, or a collection of nearby sources, using the 
higher-resolution 5-GHz data obtained with the VLA. 
Due to the complex sidelobe structure of the synthesized 
beam, nearby sources must be removed simultaneously 
from the intcrferometric data; we therefore fit all sources 
within 45" of the primary source location, roughly twice 
the synthesized beam width of the long baseline maps. 

Once we have identified all sources near the identified 
map peak which are to be removed from the data, as 
well as their morphology (compact/extended), we solve 



for source properties by fitting to the multi-pointing vis- 
ibility data. For computational expediency, we describe 
the sources as functions with analytic Fourier transforms. 
Compact sources are treated as delta functions charac- 
terized by a location, total intensity, and a spectral in- 
dex across our 8 GHz bandwidth. Extended sources are 
treated as elliptical Gaussians, characterized by a loca- 
tion, integrated intensity, spectral index across our band, 
ellipse eccentricity and angle of rotation. Parameters for 
unresolved sources are fit only to the long-baseline data, 
while those of extended sources are fit to all data. The 
best-fit models are removed from the Fourier data, and 
the mosaics are regenerated. This process is repeated it- 
eratively until there are no sources brighter than 5cr in 
the significance maps. 

In one case, the limited dynamic range of the SZA re- 
sulted in non-negligible residuals after source removal. 
For this (> 100 mjy) bright source, we remove it from 
our data, check for the greatest residual level in the vicin- 
ity of the source, and remove a section of our survey af- 
fected by the residuals which are not within our noise 
properties. This results in a hole (of 6.5' diameter) in 
our coverage where the bright source is located. 

2.6.2. Second Stage: Faint source removal 

The second stage of the fitting algorithm extracts dim- 
mer sources, relying heavily on our VLA 5 GHz follow- 
up. We develop a source catalog from follow-up data at 
5 GHz and a spatial template of these source. Begin- 
ning with the catalog, we exclude those sources already 
removed by the procedure of §2.6.1 and examine the fiux 
in the residual mosaics at the positions of all remain- 
ing 5 GHz sources. We consider any pixel whose snr is 
greater than 3 (flux > 0.5 mJy) as a source candidate. 
For these sources, we fit only for the fiux; locations are 
fixed to the coordinates from the VLA, and 5/31 GHz 
spectral indices are constrained at the lower frequency 
by the VLA fiuxes. For these dim sources, an unresolved 
source model (^-function) fit to an extended source will 
result in residual structure in the map that is within our 
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TABLE 3 
SZA Area Coverage 



Noise Value 


Area(<noise) 


(mjy) 


(dcg") 


0.22 


0.000 


0.23 


0.041 


0.24 


0.287 


0.25 


1.043 


0.26 


1.807 


0.27 


2.179 


0.28 


2.365 


0.29 


2.476 


0.30 


2.568 


0.35 


2.892 


0.40 


3.116 


0.45 


3.354 


0.50 


3.823 


0.55 


4.013 


0.60 


4.120 


0.65 


4.207 


0.70 


4.278 


0.75 


4.336 



noise properties, so we perform these fits only to the long 
baseline data and remove the corresponding flux from the 
short baseline map. 

3. RESULTS 
3.1. Survey Area 

The principal data product of our survey is the number 
of clusters detected and the area observed as a function of 
sensitivity. Equipped with this information, a means to 
translate from SZ flux to mass, a suitable mass function, 
and the number of clusters detected in our fields, we can 
estimate ag, the amplitude of mass density fluctuations 
on a scale of 8h~^ Mpc. Since the long baselines of the 
SZA are not sensitive to the scales subtended by galaxy 
clusters, the following calculations use only the source- 
subtracted, short-baseline data. 

We first calculate the sensitivity directly from our 
short-baseline sky maps (which include diffuse emission), 
and compare this to theoretical predictions. One of our 
fields (the SZA3 field) shows excess noise of 24% over 
the expectation, and there is strong evidence for a spa- 
tial correlation between the 31 GHz maps and ridges of 
dust observable in the IRAS 100/im maps. Lacking a 
suitable template for removing this foreground from our 
data, we exclude the field from our cosmological analy- 
sis. This field is the subject of a companion paper (Leitch 
et al., in prep). The rescaled, source-subtracted signifi- 
cance maps for the three remaining fields are presented 
in Figure 3. 

We calculate numerically the differential survey 
area as a function of rms noise on intervals of 
da = 0.01 mJy/beam, shown in the top panel of Figure 
4. Over most of the survey area, the noise lies between 
0.25 and 0.3 mJy/beam. In the bottom panel we present 
the integral of the top panel to indicate the total sur- 
vey area as a function of sensitivity; the same data are 
presented in tabular format in Table 3. 

3.2. Cluster Detection 

Once we have removed sources of emission and rescaled 
by the noise, we search for clusters as decrements in 
the mosaicked maps. We consider a detection any pixel 
whose amplitude is at least 5 times the rms noise level. 



Note that the Fourier-space coverage of the compact ar- 
ray was designed to match the typical cluster profile, so 
that these mosaics have already been optimally filtered 
for cluster detection. 

No clusters were identified at > 5a significance, with 
the largest decrement having a significance of 4.3cr. This 
allows us to place an upper limit on as- 

4. EXPECTED NUMBER OF CLUSTERS 

The number of clusters of mass M that we should de- 
tect, per unit noise and redshift interval, is given by 



dm dVt dz dfT 

where p{D\M,z,a) is the completeness, the probability 
of detecting a cluster of mass M at redshift z in the pres- 
ence of noise a, jj^^^f^{M, z, a^) is the mass function, 
the predicted density of clusters as a function of mass, 
redshift, and as, and ^{a) is the survey area as a func- 
tion of noise level. 

The total number of clusters we expect to detect for a 
given cosmology is then given by the integral of Equation 
2 over mass, redshift and map noise: 

Nias) = jdzjdaj -J^^ (M, .3) <iM. (3) 

Here we assume the concordance cosmological parame- 
ters from WMAP 7-year results (Larson et al. 2010). 

4.1. The Mass Function 

We assume the mass function of dark matter halos de- 
rived from cosmological simulations of fiat ACDM cos- 
mology. We adopt the fitting formula of Tinker et al. 
(2008), and assume a redshift evolution for the over- 
density parameter given by the growth factor of Viana 
& Liddle (1996). 

4.2. Completeness 

As indicated in Equation 1, the fundamental SZ ob- 
servable is the compton-y parameter. We can relate this 
observable to the cluster mass either through observa- 
tion or simulation, both of which are subject to signifi- 
cant uncertainties. SZ-effect and X-ray observations can 
be used to determine scaling relations between compton- 
y and cluster mass (e.g., Bonamente et al. 2008), but 
these observations typically comprise a small number of 
massive clusters spanning a limited redshift range, and 
often relate quantities determined at overdensity radii 
not directly comparable to those sampled by a particu- 
lar experiment. By contrast, simulations can be accu- 
rately compared to experimental details, but the corre- 
spondence between the SZ observable and cluster mass 
is highly dependent on the accuracy of the simulated gas 
thermodynamics (Kravtsov et al. 2005). 

We use the simulation of Shaw et al. (2009) (hereafter 
S09) to translate between compton-y and cluster mass. 
This simulation combines an N-body 'lightconc' simula- 
tion with a semi-analytic model for the cluster gas, with 
a significant amount of heating from feedback processes 
(i.e., AGN feedback, star formation). Gas parameters 
in S09 have also been adjusted to match X-ray observa- 
tions of low-redshift clusters. We select clusters from 40 
compton-j/ maps, each 5 degrees on a side, with M200 
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Fig. 3. — Significance (snr) maps for the SZA4, DLS, and NOAO NDWFS fields after source removal and rescaling by the noise. Note 
the region excluded in the top map due to the limited dynamic range of the SZA. We sec no dark regions of significance greater than 4.3. 



masses ranging from 2.5 x 10^^ Mq to 3.0 x 10^^ M0, 
where M200 is the mass enclosed within a radius cor- 
responding to an overdensity of 200 times the mean den- 
sity of the universe; it is this mass that we refer to in aU 
subsequent discussion. 

To calculate the probability that the SZA would detect 
a cluster of a given mass, mock observations of the y- 
maps are performed in survey mode, in a multi-pointing 
mosaic scheme. For computational expediency, we do 
not simulate observations of large areas of simulated sky, 
but select individual clusters from the catalog and cal- 
culate what fraction of them would have been detected. 
Each simulated cluster is placed at the center of a small 
field, 0.47° on a side. On this image we overlay a grid of 
pointings in a 2-3-2 hex pattern, refiecting the spacings 



in our survey. The cluster image is weighted by the pri- 
mary beam (centered at each pointing) , and the result is 
Fourier-transformed and re-sampled onto a uw-grid that 
reflects the actual coverage of the survey data. To the 
visibilities in each pointing we add Gaussian noise, with 
weights chosen such that the resulting noise in the image 
plane mosaic is uniform in a region within a 4' radius 
of the cluster. To quantify the cluster detectability as a 
function of noise, we generate maps with rms noise be- 
tween 0.2 — 0.65 mJy/beam, in steps of 0.05 mJy/beam. 

To simulate the effect of compact source contamina- 
tion, we add unresolved sources brighter than 0.05 mJy 
to the data according to the MIO distribution for regions 
further than 0'.5 from the cluster center. To account for 
cluster-source correlations we add sources brighter than 
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Fig. 4. — Area as a function of Sensitivity for the SZA survey. In 
the top panel we show the differential area as function of sensitivity 
( da = 0.01 mJy/beam). Notice the high bump near a sensitivity 
of 0.25 mJy/beam that encompasses most of our survey area. The 
smaller bump near a coverage of 0.47 mJy/beam is a result of the 
shallow coverage on half of the SZA4 field. In the bottom panel we 
present the integrated area of the survey as a function of sensitivity. 



0.01 mJy according to the distribution from Coble et al. 
(2007) for regions near the cluster center. To mimic our 
source-extraction algorithm, we generate a 5 GHz flux 
for each source from the spectral index distribution of 
MIO, and create a 5-GHz source catalog of all sources 
brighter than 0.35 mJy (the detection threshold for our 
VLA 5-GHz source catalog). Sources are extracted from 
the simulated maps exactly as described in §2.6, with the 
exception that we remove a perfect model from the sim- 
ulated data rather than a fitted fiux and spectral index, 
in the interest of computational speed. We have verified 
that this does not systematically bias our results by com- 
paring the fitted flux of sources in these simulations with 
their original values, i.e., the presence of noise introduces 
a random error into the fit fiux of the source, yet over 
the many realizations performed we recover the true fiux 
of the source. 

Once the identified sources are removed from the mock 
observation, we search for the peak decrement in the 
short-baseline mosaics, within the region of uniform noise 
in our mock observation (namely 4'). We identify clus- 
ters in the simulated data as described in §3.2. 

To capture the inherent fiux distribution for clusters 
of similar mass, as well as the redshift evolution of clus- 
ters, clusters were observed in nine different mass ranges 
(M,„i„ = 1.5 X 10^4 ^ iQii ]^^) gjgj^^ 

shift bins (zmin = 0.1, Az = 0.2). The number of clus- 
ters observed in each mass range and redshift bin varied 
from 2 to 200; some bins contained no clusters, namely 
the high-mass, high-redshift bins. We further select only 
clusters which are not within 0.23° of the edge of the 
simulated maps, and which are not within 4 arcminutes 
of clusters of greater or comparable masses. The first 
cut ensures we do not introduce artifacts associated with 
field edges, and the second ensures that our observations 
are of clusters whose masses are well-defined, i.e., whose 
simulated masses are not corrupted by the presence of 
secondary clusters. For each mass and redshift bin, we 
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Fig. 5. — The mass limit of the SZA survey calculating us- 
ing the area as a function of noise, and cluster detectability from 
simulations of S09 (square points) and W03 (See §6). The points 
correspond to the M200 mass of a cluster which would have been 
detected at least 50% of the time, and are offset from each other 
in redshift for ease of presentation. The differences between values 
are a result of the W03 simulations not including AGN feedback 
or "pre-heating" in calculating the compton-j/ parameter. 

iterate over 200 distinct realizations of clusters, noise and 
compact sources. 

To calculate the resulting completeness, we determine 
the fraction of time the clusters in a given bin were de- 
tected at significance greater than 5a. 

4.3. Survey Mass Limit 

Equipped with the survey completeness (§4.2) and the 
area as a function of noise (§3.1), we can calculate the 
area as a function of mass and redshift. As the SZ effect 
depends most directly on cluster mass, this can be used 
to calculate the redshift-dependent mass limit of our sur- 
vey, which we present in Figure 5 for the S09 simulations. 
For each redshift bin, we present the M200 mass above 
which our survey is 50% complete. As expected, the mass 
limit decreases towards higher redshifts, as clusters be- 
come more compact and therefore more easily detectable 
by the SZA. 

In Figure 5, we also present the limiting mass cal- 
culated from a different simulation. White (2003), to 
demonstrate the effect of gas physics in our complete- 
ness calculations. The differences between the simula- 
tions and their effect on our results are discussed in detail 
in §6. 

5. CONSTRAINT ON 0-8 

To calculate the value of cts which is most consistent 
with the SZA survey, we address the question: what is 
the probability of a erg value given the number of ob- 
served clusters, P{(Js\N)7 A simple invocation of Bayes' 
theorem yields: 

Pias\N) cx P{N\as)P{<7s) (4) 

where P{N\as) is the probability of N detections given 
a value of cts and Piers) is the prior on trg, which we take 
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Fig. 6. — Expected number of clusters (and the probability of 
non-detection of clusters) in the 4.4 square degree SZA survey, as 
a function of as, using WMAP cosmology for all other parame- 
ters and cluster detectability from S09 simulations. The solid line 
represents the expected number of clusters, with the corresponding 
scale presented on the right axis. The dashed line is the probability 
of detecting zero clusters, with the corresponding scale on the left 
axis. 



to be uniform. 

The number of clusters detected is a Poisson process, 
with a well-defined distribution given by 



P{k\X) 



-Aye 

kl 



(5) 



where P{k\X) is the probability that a process with ex- 
pectation A occurs k times. As we detect no cluster can- 
didates in the SZA survey (at a significance greater than 
5a), we have k=0, whence 



P{as\N) = e 



(6) 



where N^a^) is the result from Equation 3, shown pic- 
torially in Figure 6. Note that in determining N(as), 
we adopt concordance cosmological parameters from 
WMAP 7-year results (Larson et al. 2010), and calculate 
for completeness determined from the S09 simulations. 
Figure 6 (via Equation 6) can be used to compare the 
relative likelihood of different values of erg. For example, 
we expect a non-detection in the SZA survey 5% of the 
time if erg = 0.97. 

We integrate Equation 6 to generate the probability, 
given our data, that erg exceeds a given value; the result 
is presented in Figure 7. We see that for the S09 simula- 
tion, as values of greater than 0.84 are ruled out at 95% 
confidence. 

In calculating this limit, we have made the conven- 
tional assumption of a "non- informative" prior (i.e., uni- 
form for (Tg > 0). While we consider this approach 
conservative, in that it captures what the; SZA data can 
say about the value of cxg under a minimal set of prior 
assumptions, note that it will also lead to the tightest 
possible (for uniform priors) constraint on erg. If for ex- 
ample, the prior were truncated to cg > 0.5, our 95% 
limit would shift to erg < 0.87. Sensitivity to the prior is 




Fig. 7. — Plot of the probability of that as exceeds a given value 
for completeness derived using the S09 simulations, with the 1-, 2-, 
and 3-<T confidence limits presented. 

typical for data sets with unequal power to discriminate 
among allowable hypotheses; in this case, the small area 
of the SZA survey can strongly rule out high values of 
as, but provides little or no power to discriminate against 
low values of erg, to which our prior nonetheless assigns 
equal weight. 

6. SYSTEM ATICS 

The conversion from cluster mass to SZ observable is 
potentially the largest systematic uncertainty in cosmo- 
logical analyses of SZ surveys. Because constraints on 
(Tg are determined by comparison with simulations, they 
are necessarily sensitive to the assumptions that under- 
lie each simulation's modeling of the ICM. As a rough 
estimate of the importance of model assumptions to our 
limit on erg , we repeat our analysis using the simulation of 
White (2003) (hereafter W03), a high- resolution N-body 
simulation following the semi-analytic method of Schulz 
k White (2003). The salient difference between the S09 
and W03 simulations is that S09 incorporate AGN and 
supernovae feedback while W03 ignore cold gas and star 
formation. As a result, clusters in the W03 simulation 
tend to be more compact (and therefore more detectable 
by the SZA), with higher central SZ-decrements than 
equivalent clusters (by mass) in S09. This leads to a 
higher completeness for W03 clusters of a given mass, a 
lower mass limit (See Figure 5), and a stronger constraint 
on erg (a 95% confidence limit that is lower by 0.04). We 
choose to present the limit derived from S09 in §5 not 
only because it is the more conservative of the two, but 
also because the S09 simulation represents a more real- 
istic model of ICM physics; we nevertheless caution that 
significant uncertainties in the modeling of cluster gas 
physics remain. 

In the completeness calculation described in §4.2, we 
model the correlation between clusters and radio sources 
by using the Coble et al. (2007) source distribution for 
the inner regions of clusters, and the MIO distribution 
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for field sources. However tfie Coble et al. (2007) dis- 
tribution was determined from observations of the most 
massive clusters, and can bias our completeness low if 
sources are over-represented in these objects relative to 
the lower-mass clusters to which the SZA is sensitive. 
To bracket the magnitude of this eflFect, we repeat the 
completeness calculations described in §4.2 without us- 
ing the Coble et al. (2007) over-densities, i.e., using the 
source counts obtained only from field sources presented 
in MIO. Neglecting the higher density of sources towards 
galaxy clusters increases our completeness on the order 
of 12%, which due to the steepness of the dependence of 
cluster counts on (jg would lower our limit by 0.04. 

In §2.6 we mentioned that for one bright source, the 
limited dynamic range of the SZA introduces a hole in 
our coverage at the source location. If the correlation 
between clusters and radio sources had not properly been 
taken into account in calculating our completeness, this 
would bias our completeness low if bright sources are 
preferentially associated with galaxy clusters. However, 
we simulate both the correlation and the corresponding 
reduction in completeness entailed by missing clusters 
associated with sources that we cannot subtract; cutting 
such sources out of our maps therefore corresponds to a 
simple reduction in the survey area, and should have no 
effect on our limit. 

Although we see no evidence of soiirces which are re- 
solved by the long-baseline data, as discussed in §2.6, we 
fit the 12% of sources deemed to be extended at 5 GHz 
to a combination of the long and short-baseline data 
(whereas sources deemed a priori to be unresolved were 
fit only to the long-baseline maps) . If any of these sources 
happened to be associated with a cluster, the cluster 
decrement present in the short-baseline data would re- 
duce the fitted flux of the source, resulting in residual 
source flux in the map that would bias against detect- 
ing the cluster. To quantify this effect, we re-calculate 
the completeness with a slight modification in 12% of 
the realizations, namely that we do not remove the true 
flux of the source, but instead remove the map flux at 
the source location in the combination of short and long 
baseline maps. As the majority of sources are not as- 
sociated with clusters, this test can be thought of an 
upper limit to any potential reduction in our complete- 
ness. Properly accounting for this effect would raise our 
limit on ag by at most 0.01. 

An additional source of uncertainty unrelated to our 
analysis methods is the clustering of large-scale struc- 
ture. We know that galaxy clusters are not evc;nly dis- 
tributed throughout space, but form preferentially along 
fllaments. As a result, the number of clusters seen in 
small fields such as those observed by the SZA do not 
follow a Poisson distribution. In particular, since most 
lines of sight will sample the voids between filaments, this 
leads to an increased probability of detecting fewer mas- 
sive halos than the Poisson average. To quantify the im- 
pact on our results, we begin with the 100 square degree 
W03 simulations. We select regions similar in shape, size 
and noise properties to each of our fields from the sim- 
ulated maps, and fold in the calculated completeness to 
estimate the number of clusters we should detect. We do 
this for > 10'* realizations of each of the three SZA fields, 
varying the location and orientation of the field over the 
sky-map, to generate a distribution of expected numbers 



of clusters, and compare this to the Poisson prediction 
for the equivalent field size and input erg- The result is a 
10% increase in the probability of a non-detection for an 
input (Tg of 0.9 over the Poisson prediction, leading to an 
underestimate of the true erg. Accounting for this effect, 
assuming similar clustering over a range of as, would 
raise our limit by at most 0.04. 

Lastly, as described in MIO, the calibration of the SZA 
data is tied to the modeled fiux of Mars, which we esti- 
mate is accurate to < 10%. A shift in the fiux of Mars 
corresponds to a simple rescaling of the map noise, and 
consequently a shift in the completeness. As noted above, 
a change in the completeness of < 10% corresponds to a 
shift in our limit on as (in either direction) of < 0.04. 

7. DISCUSSION 

From 2005 to 2007, the SZA performed a 6.1 square 
degree survey for clusters of galaxies via their SZ effect 
at 31 GHz. In one of the fields there is evidence for large- 
scale dust-correlated emission; this field was excluded 
from the analysis presented here. Of the remaining 4.4 
square degrees of the survey suitable for cosmological 
analysis, we estimate that the survey is 50% complete to 
a mass of M200P 6 X 10" Mq, averaged over redshift. 
By comparison with simulations, we place an upper limit 
on the value of as of 0.84 (+0.07) at 95% confidence, 
where the uncertainty refiects calibration and systematic 
uncertainties discussed in §6, excluding the error associ- 
ated with the simulated cluster gas physics. Although 
this last uncertainty is potentially the dominant one, to 
properly quantify it requires a calculation of the com- 
pleteness for a wide range of simulations with realistic 
gas models, which is beyond the scope of this paper. 

Our limit on erg is consistent with recent results from 
SZ surveys performed over larger areas of sky, such as 
with the South Pole Telescope (Vanderlinde et al. 2010), 
and with determinations of erg from gravitational lens- 
ing and X-ray cluster surveys (Smith et al. 2003; Allen 
et al. 2003). In addition, it is consistent with determina- 
tions of as from CMB anisotropy measurements, namely 
those of WMAP (Dunkley et al. 2009), the South Pole 
Telescope (Lueker et al. 2010), and the Atacama Cos- 
mology Telescope (Fowler et al. 2010). Our constraint is 
also in agreement with that of Sharp et al. (2010), based 
on CMB anisotropy measurements with the SZA itself. 
Although the data were collected with the same instru- 
ment, we note that both the data sets and the analyses 
of this paper and of Sharp et al. (2010) are completely 
independent. 
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